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C/3 

O . Abstract: The fluctuations in the elastic light scattering spectra of normal 

and dysplastic human cervical tissues analyzed through wavelet transform 
based techniques reveal clear signatures of self-similar behavior in the 
spectral fluctuations. Significant differences in the power law behavior 
ascertained through the scaling exponent was observed in these tissues. The 
strong dependence of the elastic light scattering on the size distribution of 

jl^.,^, the scatterers manifests in the angular variation of the scaling exponent. 

(~| ' Interestingly, the spectral fluctuations in both these tissues showed multi- 

0^, fractality (non-stationarity in fluctuations), the degree of multi-fractality 

being marginally higher in the case of dysplastic tissues. These findings 
using the multi-resolution analysis capability of the discrete wavelet trans- 

^ ' form can contribute to the recent surge in the exploration for non-invasive 

0^ , optical tools for pre-cancer detection. 
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1. Introduction 

The use of optical techniques for the study of biomedical systems is a rapidly developing field 
that has seen a dramatic expansion in the recent years, partly due to tremendous progress in the 
field of lasers, fiber optics and associated technologies. Both medicine and biotechnology re- 
quire appropriate instrumentation to analyze and monitor biological systems for deviations from 
normality. Optical methods, due to their non-invasive nature, are providing novel approaches for 
medical imaging, diagnosis and therapy. Considerable efforts have been made in the recent past 
to use optical techniques such as fluorescence spectroscopy [1, 2, 3, 4], Raman spectroscopy [5] 
and elastic scattering spectroscopy [6] for quantitative and early diagnosis of various diseases. 
Several optical imaging techniques like coherence gated imaging, polarization gated imaging 
and diffuse optical tomography are also being actively pursued for obtaining high resolution 
(micron scale) images of biological objects and their underlying structure [7, 8, 9, 10, 11]. 

For optical diagnosis, elastic and inelastic light scattering spectra from tissues are exploited. 
The in-elastically scattered light (via processes like fluorescence and Raman) contain useful 
biochemical information about the sample that can be employed for probing subtle biochemical 
changes as signatures of disease progression. On the other hand, elastically scattered light from 
biological tissues contain rich morphological and functional information of potential biomed- 
ical importance [6, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]. Both the angular 
and wavelength dependence of the elastically scattered light from tissue can be analyzed to 
extract and quantify subtle morphological changes taking place during progression of a disease 
[14, 15, 16, 17, 18, 21, 22], and thus may be explored as a sensitive tool for early diagnosis. 
This would however involve appropriate modeling of light scattering in complex random media 



like tissues, and the development of suitable approaches to extract/interpret the morphological 
information contained in the elastic light scattering signal. 

The spatial fluctuation of the refractive index in biological tissues arising from scatterers 
ranging in sizes from a few nanometers to a few micrometers give rise to elastic scattering 
[6, 20]. The lack of sufficient knowledge about the complex dielectric fluctuations in the tissues 
pose a formidable problem in the exact modeling of light scattering. Nevertheless, several ef- 
forts have been made in the recent years using electromagnetic (EM) theory based approaches 
like Mie theory and Born approximation to model and understand the scattering process from 
biological tissues [26, 27, 28, 29, 30]. It has also been shown that the refractive index fluctua- 
tions in biological tissues are fractal in nature which can be used to understand the structural 
changes in tissues induced by diseases [16, 27, 31, 32, 33]. 

Since the tissue morphology dependent refractive index fluctuations are recorded in the elas- 
tic scattering spectra [34], analysis of elastically scattered spectral fluctuations using sophisti- 
cated fluctuation analysis tools might facilitate extraction and quantification of subtle morpho- 
logical changes associated with early stages of cancer The scaling behavior which is generally 
assumed to be global (mono-fractal), has been shown to manifest in the local fluctuations in 
various physical processes [35, 36] and has been characterized using Multi-Fractal De-trended 
Fluctuation Analysis(for example see [37]). Wavelet Based Multi-Fractal Detrended Fluctua- 
tion Analysis (WB-MFDFA) is one other state-of-the-art technique that can be used for extract- 
ing and quantifying the self similarity at varying length scales associated with the structural 
changes associated with cancer progression due to the inherent use of fractal like transforma- 
tion kernels. 

Wavelet transform due to it's multi-resolution analysis capability using the Daubechies' ba- 
sis which extract the polynomial trends (for example, Db-4 and Db-6 extract the linear and 
quadratic trends respectively) has been shown to characterize the scaling behavior and self- 
similarity of empirical data sets quite faithfully [38, 39]. Indeed, it has been initially explored 
to analyze tissue fluorescence spectra in an attempt to distinguish between normal and dysplas- 
tic tissue [40, 41, 42, 43, 44]. In this work, we employ this multi-resolution property of wavelets 
to ascertain the changes in the self-similarity of dysplastic human cervical tissues as opposed 
to healthy human cervical tissues by analyzing the esoteric nature of the fluctuations in tissue 
light scattering spectra. 

This article is organized as follows: Fourier and power spectrum analysis is reviewed in 
2.1, discrete wavelet transform in 2.2, wavelet based power law analysis in 2.3, wavelet based 
multi-fractal de-trended fluctuation analysis in 2.4 and correlation based analysis in 2.5. Sec. 
3 describes the experimental methods for light scattering measurements from tissues. Sec. 4 
deals with our findings from the analysis and contains a discussion of the same in the context 
of the differences between the normal and dysplastic samples. In Sec. 5 we conclude with 
the prospect of pre-cancer detection using light scattering techniques combined with novel 
fluctuation analysis methods. 

2. Theory 

2.1. Fourier analysis and power law spectrum 

Fourier Analysis has traditionally been a preferred tool for analysis of experimental data sets. 

Here, we just briefly review the Discrete Fourier Transform (DFT) and its power spectrum. For 

a data set x(n),n = {1,2,. . .}, the DFT is a linear transformation over an orthogonal basis given 

by: 

^-1 / 2ki \ 

x(^)= ^^(«)exp -— foz , ^e[0,A^-l] (1) 

n=0 V ^ / 



and it's power spectrum 
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n=Q V 



(2) 



is known to follow a power law behavior for self-similar processes 

P{k)^k". (3) 

The power law coefficient a is related to the Hurst scaling exponent i/ G (0, 1), and the frac- 
tal dimension Df which identify the nature of fractal behavior by a = 2H +1=3— Df [35]. 
Though this analysis assumes a global scaling behavior, due to the presence of inhomogeneities 
in the intra-cellular structure, the scaling behavior turns out to be locaUzed in the spatial fre- 
quency domain. In multi-fractal analysis, a local Hurst exponent is calculated which quantifies 
the local singular behavior and provides useful and minute information which is hidden be- 
tween different scales. We therefore applied a more general wavelet based approach for the 
multi-fractal analysis. 

2.2. A brief review of wavelet transforms 

Wavelet Transforms [45, 46, 47, 48], in both the forms of Discrete and Continuous transforms 
in the recent years have emerged as an invaluable tool in the field of data analysis and interpre- 
tation. Here, we will briefly review the Discrete Wavelet Transform (DWT). 

In DWT, two functions, namely, (j){n) and </(«), called the father and the mother wavelets 
respectively, form the kernels. They satisfy the admissibility conditions: f(p{n)dn < °o, 
J\(/{n)dn = 0, J(j)*{n)\i/{n)dn = 0, J\(j){n)\^dn — J \\j/{n)\^dn = 1. The scaled and translated 
versions of the mother wavelet \(f{n) are called the daughter wavelets 

V/s,,„(n) == 2'/V(2'n - m), m G M, s e Z+ , (4) 

which form a complete set and differ from the former in terms of their height and width. At 
the 5''' scale, the height and width of the daughter wavelet are 2* and 2^' ^ of that of the mother 
wavelet respectively, s and m are the scaling and translation parameters. 
The DWT for a function f{t) is given by, 

I 

f{t)= Y^ a,n<p,n{t)+ 52 Y^d,,„,\if,jn{t), (5) 

m=— oo m=— ooi>0 

the coefficient am (called the approximation coefficient) extracts the trend and c/^ ,„ (called the 
detail coefficient) extracts the details or fluctuations from the signal. Here, / = [log (A^)/ log (2) J 
is the upper bound for taking the maximum scale for analysis above which the edge effects 
corrupt the wavelet coefficients giving rise to spurious results. This mathematical artifact is 
explained by the cone of influence [47]. 

2.3. Wavelet Fluctuation Power Law Analysis 

The Daubechies' family of wavelets are made to satisfy vanishing moments conditions, and 
hence are able to isolate polynomial trends from a time series (see [45]). The different wavelets 
in this family isolate trends of different polynomial orders, for example, Db — 4 isolates a mono- 
mial trend while Db — 6 isolates a quadratic trend. We first obtain the profile from the fluctu- 
ations by taking their cumulative sum: Y{yi) = ^"j^iX(m),n G {1,2,...,A^— 1}, where, Y and 
X are profile and fluctuation signal respectively and A^ is the data length. In this method, we 
first obtain the fluctuations at every scale by a wavelet reconstruction taking the approxima- 
tion coefficients (using Db-4 wavelet) and subsequent subtraction from the signal. A flowchart 




Fig. 1. Schematic representation of the fluctuation extraction algorithm (adapted from [49]). 



for this fluctuation extraction is shown in Fig. 1. We then fit the Fourier power spectrum of 
these fluctuations to obtain the power law exponent as a function of scale, i.e. a = a{s) for the 
identification of the short and long term correlations. This technique sheds light on the scaling 
behavior of the fluctuations in different spatial frequency regimes. 

2.4. Wavelet Based Multi Fractal De-trended Fluctuation Analysis 

The Multi Fractal De-trended Fluctuation Analysis (WB-MFDFA) algorithm proposed by Man- 
imaran et al. [38, 39] has been employed gainfully to extract the multi-fractal nature of a variety 
of physical systems. In order to apply the WB-MFDFA algorithm, we then obtain the fluctu- 
ations using the extraction algorithm shown in Fig. 1. The asymmetric nature of the wavelet 
function and the edge-effects (due to the cone of influence) encountered during the convolu- 
tion can affect the precision of the fluctuations. Hence, we repeat this procedure on a reversed 
profile and then take the average to get the fluctuations at every scale which are denoted by F^. 
We used s to represent the wavelet scale so as not to confuse with the segment length s which 
is related to the wavelet scale s by the number of filter coefficients for a given wavelet. These 
fluctuations are then segmented into M^ non-overlapping sections such that Ms — \N /s\ , where, 
s and A^ are the window size and the length of the fluctuations respectively. Subsequently, we 
obtain the q"^ order fluctuation function Fq{s) by 



F,{s) = 



2M, 



ir I {^'(-'^)} 



■' m=l 



q^Q 
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where the order of moments q can and both positive and negative integers. The fluctuation func- 
tion for self similar processes follows a scaling law, the scaling function given by /^y(i) ~ s''^'^' . 
We should note here that the smaller fluctuations in the light scattering spectrum will be influ- 
enced by the negative q values, whereas the positive values of q will impact the larger fluctu- 
ations. The scaling function li{q) calculated at ^ = 2 corresponds to the Hurst exponent [35]. 



// = 0.5 represents uncorrelated (white noise, /") or brown noise /^^, while H >0.5 represents 
long range correlations or persistent behavior. H < 0.5 reveals short range correlations or anti- 
persistent behavior Mono-fractals are scale independent and hence their h{q) is independent 
of q. They can be characterized by a single parameter like the fractal dimension. However, for 
multi-fractals, the h{q) is not independent of q and they require a more complex function like 
the singularity spectrum for its characterization [36]. 

The classical multi-fractal scaling exponent z{q) defined by standard partition function based 
formalism [50, 51] is related to the Hurst exponent h(q) by T{q) = qh{q) — 1. The multi- 
fractality can also be characterized by using the singularity spectrum /()3), which is related 
to z{q) through a Legendre transform; j3 = d/dq[T{q)] and /(j3) ^ qj5 — t(^); where j3 is the 
singularity strength or Holder Exponent [37]. /(j3) denotes the dimension of the subset of the 
series to be characterized. j3, /(j3) and h{q) are related by 

f{P)^q[l5-h{q)] + l. (7) 

/3 will have a constant value for mono-fractal series, whereas for multi-fractal series j3 values 
will have a distribution. Just like in the case of h{q), a constant t(^) will signify a mono- 
fractal, while a multi-fractal will depend on the order of the moments. The singularity spectrum 
is broader when correlations are present in the series and if the correlations are absent or lost, 
then the singularity spectrum becomes narrower 

2.5. Correlation based analysis 

The correlation matrix analysis method has been extensively used to study various physical and 
biological systems (for example see [52]). For two variables Xj E X and yj G Y, the correlation 
matrix is given by: 

where, (. . .) is the mean. It is well known that multi-fractality is related to the domain struc- 
ture formation in the correlation matrices and we present the correlation matrix analysis in the 
wavelength domain to support the wavelet based multi-fractal de-trended fluctuation analysis. 

3. Experimental materials and methods 

Pathologically graded (CIN-I or dysplastic) unstained sections of human cervical tissue and 
their normal counterparts were frozen and vertically sectioned (thickness ~ 5jUm, lateral di- 
mension ^ Amm X 6nim) on glass slides were used for light scattering measurement and analy- 
sis. A schematic representation of the experimental setup is shown in Fig. 2. The angle resolved 
elastic scattering spectra were recorded using a goniometric arrangement for the spectral range 
X : AQQnm - 800«m and angular range 0:10°- 150° at 10° intervals. White light output from 
a Xe-lamp (Newport USA, 50-5000 W) was collimated using a combination of lenses and were 
made incident on the sample kept at the center of the goniometer. The spot-size incident on the 
sample was controlled by a variable aperture {^ \mm). The scattered light from the sample was 
collimated by a pair of lenses and was then focused into a collecting fiber probe, the distal end 
of which was coupled to a spectrometer (Ocean Optics HR2000). The recorded spectra at each 
scattering angle was normalized by the lamp spectra. 

4. Results and Discussions 

Typical elastic light scattering spectra recorded from normal and dysplastic tissues are shown 
in Fig. 3. Representative spectra for forward and back-scattering are shown for Q = 10° and 




Fig. 2. (Color Online) Schematic representation of the experimental setup for the light 
scattering measurement. The sample mounted on a goniometer is illuminated by a white 
light coUimated from a Xe-source and the scattered light is then collimated using two lenses 
(L2 and L3) and recorded using a Fiber-Optic Spectrometer. The Of represents the forward 
scattering angle, while 0/, is the backward scattering angle. The forward and backward 
scattering process are represented by red and blue lines respectively while the Fiber-Optic 
Spectrometer setup is shown in green. < < 180° is the scattering angle at which the 
spectra were recorded. We must note in the gray shaded region extending from 80° to 1 10°, 
the scattering spectra was not recorded due to geometry induced artifacts. 



= 150° in Figs. 3(a) and (b) respectively. The broader structures are possibly due to larger 
sized scatterers, whereas the smaller features hidden in the spectra are the fluctuations due to 
the smaller sized scatterers [13]. We can see that the spectrum flattens as we move towards the 
backward scattering angles as compared to the forward scattering angles, which signifies the in- 
volvement of larger sized scatterers, such as the nucleus, in forward scattering and smaller sized 
intracellular organelles in the backward scattering. Such angular dependence of elastic scatte- 
ring spectra from tissues have been extensively studied both experimentally and theoretically. 
These studies have indeed shown that the contribution of larger scatterers (like cells, nuclei) 
are more dominant in the spectra recorded at forward angles, whereas the spectra recorded at 
back- scattering angles are typically more influenced by the smaller sized scatterers [24]. 

Although it is difficult to make a one-to-one correspondence between the intensity fluctua- 
tions and local refractive index fluctuations in the morphological structures, some information 
about the morphology can still be inferred from the light scattering spectral fluctuations. In- 
deed, it has been shown that with Born approximation, the light scattering spectrum can be 
represented as a Fourier transform power spectrum of the local refractive index fluctuations 
(see Sec. 2. 1 for details). Thus, the observed spectral fluctuations contain information about the 
corresponding refractive index fluctuations in the Fourier domain as discussed below. 

The results of the Fourier analysis on the light scattering spectra for different scattering 
angles is shown in Fig. 4. In the forward scattering region (10° — 50°), we observe that the 
dysplastic samples show ~ 0.7 < a < 1.2, in the same region, however, the normal tissues 
show 0.4 < a < 0.8 (where a is the power law coefficient). The higher values of the power- 
law coefficient a corresponds to higher values of the Hurst scaling exponent H, indicative of 
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(a) Data at 10°. 
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Fig. 3. Recorded elastic light scattering spectra at (a) 10° and (b) 150° from both normal 
and dysplastic tissues. 



"coarseness" of the fluctuations tending towards sub-fractal behavior; whereas at higher angles 
(1 10° — 150°), a values are observed to be smaller in dysplastic samples than that for normal 
samples in the same region indicating more "roughness" in the fluctuations, signifying a trend 
towards extreme fractality. 

The higher a values in the angular range (10° < < 70°) for the dysplastic tissues as shown 
in Fig. 4, due to the more dominant contribution of the large scale scattering structures. This 
possibly can be related to the proliferative nuclear morphology in dysplasia [13]. This fol- 
lows because such morphological changes are expected to lead to coarser fluctuations in light 
scattering spectra as is evident from the higher value of H in the dysplastic tissues in the for- 
ward scattering range. On the other hand, the lower values of a in the back-scattering angles 
(110° < < 150°) indicates the predominance of small scale heterogeneity in the dysplastic 
tissues. It is also pertinent to note that around 60°, there is no distinct difference, possibly due 
to the simultaneous contribution of different sized scatterers in this angular range resulting in 
"lumped" effects. 
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Fig. 4. Variation of the scaling exponent a as a function of the scattering angle as obtained 
from the Fourier analysis. 



In Fig. 5, we present the plots of power-law scaling exponent a as a function of the scale 



for normal (Fig. 5(a)) and dysplasia (Fig. 5(b)) for forward scattering angles 40° and 50° and 
for backward scattering angles 140° and 150°. This analysis was performed following the dis- 
cussion in Sec. 2.3. It must be noted that the difference in the absolute value of a for Fourier and 
wavelet analysis arises from the fact that, in Fourier analysis, we analyze the signal itself while 
in wavelet analysis, we obtain the a values from the fluctuations which can be thought of as 
higher derivatives of the signal. We observe that the a values show a strong dependence on the 
scattering angle 0. In both Fig. 5(a) and Fig. 5(b), we observe the presence of a broad spectrum 
of processes such that 1 .4 < a < 3.2. This arises possibly from the contribution of differently 
sized scatterers as has been mentioned earlier We must however note that the behavior of the 
normal and dysplastic tissues in power law regime turn out to be very different. For example, 
in the forward scattering region, at 50°, we observe that the normal tissues show an average 
a = 2.6, while the dysplastic tissues show an average a = 2.0 over all scales. Similarly, in the 
backward scattering region, we obtain a = 1.8 and a = 3.1 for normal and dysplastic tissues 
respectively at 150°. In the other regions also, we find stark differences between normal and 
dysplastic samples in the a values averaged over all scales. 
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Fig. 5. Variation of the power law exponent a with the wavelet scale s at representative 
forward and backward scattering angles for (a) normal and (b) dysplastic tissues. 

As discussed in Sec. 2.4, we have depicted the Hurst exponent as a function of the scattering 
angle 6 in Fig. 6; where we notice that though for both the normal and dysplastic tissues show 
a "near-random" behavior corresponding to // = 0.5 (// = 0.5 is indicative of Brownian behav- 
ior), we observe that the dysplastic samples show a marked departure from the normal samples 
in their trend. As observed for a{6) in Fig. 4, H has a higher value in the forward scattering 
range 10° < < 70° while the trend reverses in the backward scattering range 120° < < 150°. 

The singularity spectrum /(j3 ) is a quantitative indicator of the exact nature of the self- 
similarity and its width represents the strength of the multi-fractality. In Fig. 8, we have shown 
the /(j3) as obtained from Eq. 7 as discussed in Sec. 2.4. We observe that the dysplastic samples 
have a higher multi-fractality than the normal samples, indicated by the width of singularity 
spectrum. It must be noted that for a mono-fractal, the singularity spectrum is similar to a 
Gaussian with a very small variance. This is consistent with our observations of Fig. 7, where, 
the plots of h{q) vs q for normal and dysplastic tissues at a few representative forward and 
backward scattering angles are shown. For mono-fractals, the h{q) is independent ofq. In Fig. 7, 
we observe that for the normal samples, at angles 50° — 70°, the dependence of h{q) on q is very 
small {Ali{q) ^ 0.2), which implies a trend towards mono-fractality. The normal samples thus 
show a wide range of variation from mono-fractals to multi-fractals. However, for dysplasia. 
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Fig. 6. The extracted Hurst parameter, H = h(q = 2) as a function of the scattering angle 6 
for normal and dysplastic tissues. 
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Fig. 7. The scaling function h{q) at different forward and backward scattering angles for 
(a) normal and (b) dysplasia. The weaker q dependence of h{q) for normal samples in the 
forward scattering angles (50° — 70°) is indicative of a possible mono-fractal trend, while 
the stronger dependence of the scaling function on the order of moments for dysplastic 
samples is indicative of a multi-fractal trend. 



the h{q) dependence on q is high (Ah{q) ^ 0.7 ) for all scattering angles, which indicates that 
the dysplastic samples show more multi-fractality than the normal samples. 

The spectral correlation matrices are shown in Fig. 9 (following Sec. 2.5). It is clear from 
Fig. 9(a) and 9(b) that though the normal tissues do not show any distinct correlation sectors in 
this domain, other than the expected correlation that would occur at neighboring wavelengths; 
dysplastic tissues show the presence of three dominant sectors. It must also be noted that the 
range of correlation increases from 0.70 — 1.00 for normal to 0.97 — 1.00 for dysplasia. This 
indicates a higher correlated behavior of the dysplastic tissues than the normal tissues in addi- 
tion to domain formations in the spectral range. This possibly arises due to the fact that during 
dysplastic progression, the homogeneous cell morphology gives way to a more fragmented and 
heterogeneous structure. 
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Fig. 8. The singularity spectrum /(jS) plotted against j6 at all scattering angles 9. (a) shows 
the singularity spectrum for normal samples while (b) shows the singularity spectrum for 
dysplastic tissues. 
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Fig. 9. (Color Online) The correlation matrices in the wavelength domain for (a) normal 
and (b) dysplastic samples. 



5. Conclusion 

In conclusion, we have applied a combined Fourier based and discrete wavelet based analy- 
sis on the fluctuations extracted from the elastic light scattering spectra of normal and dys- 
plastic human cervical tissues. This approach clearly revealed otherwise hidden signatures of 
self-similarity in spectral fluctuation for both normal and dysplastic tissues, with significant 
differences in the nature of self-similarity. Fourier analysis of these fluctuations indicated the 
existence of multi-fractal nature and was further confirmed by WB-MFDFA. Dysplastic tis- 
sues showed marginally higher multi-fractality over the entire angular region compared to their 
normal counter-parts. The scaling exponent was observed to have angular dependence, possi- 
bly arising from the size distribution of the scatterers present in such complex systems. Note 
that this novel fluctuation analysis approach was initially explored on elastic scattering spectra 
recorded from a limited number of tissue samples, and the results were qualitatively similar. A 
more systematic study on statistically significant number of samples is currently underway, and 
the results will be reported in the near future. Never-the-less, early indications show promise of 
this approach for quantification of the morphological alterations associated with pre-cancers. 



The understanding gained from such analysis on the muhi-fractal nature of fluctuations in tis- 
sues may ultimately lead to the development of non-invasive optical tools for pre-cancer detec- 
tion. 



